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Abstract 

^ i We use linked cluster series expansion methods to estimate the values of 

, various short distance correlation functions in S = 1/2 Heisenberg antiferro- 

magnets at T = 0, for dimension d = 1,2,3. The method incorporates the 
\ possibility of spontaneous symmetry breaking, which is manifest in d = 2, 3. 

O \ The results are important in providing a test for approximate theories of the 



antiferromagnetic ground state. 
> ■ PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



This paper deals with the problem of calculating correlation functions, at T = 0, for the 
S —1/2 Heisenberg antiferromagnet 

i/ = 2J^S,-S,- (1) 

where the sum is over all nearest neighbor pairs. We consider explicitly the linear chain 
{d = 1), the square lattice (d = 2) and the simple cubic lattice (d = 3). The correlation 
functions (correlators) are defined as 

C(r) = 4(Si • Si+r)o = Ci{r) + 2Ct(r) (2) 

where the average is a ground state expectation value, r is the distance between sites in units 
of the lattice parameter, and the factor 4 is included for convenience. It is also convenient 
to separate the correlator into a longitudinal part C;(r) and a transverse part Cj(r) 

Ci{r) ^ 4(5767+,)o (3a) 
a(r)^2(5f5^, + 5|'5|V.>o 

= {^i^i+r + Stvr)o (3b) 

Ci and Ct will differ if the Hamiltonian is generalized to include Ising anisotropy 

H^2JJ2 [S-S] + X{S^S^ + S^S])] (4) 

(ij) 

as we shall do, or if the ground state of the isotropic Hamiltonian exhibits spontaneous 
symmetry breaking. We shall see that this occurs in dimension d > 1. 

The correlators characterize the nature of the ground state of the system, and hence 
an accurate knowledge of their values can be important for testing approximate analytic 
theories. Surprisingly, apart from the 1-d case, knowledge of their values is limited. 

We use the method of linked cluster expansions in which the Hamiltonian is written as 

H ^Ho + XV (5) 



with the Ising part taken as the unperturbed Hamiltonian and the remainder as a perturba- 
tion. To improve the convergence of the series, we also include a local staggered field term 
tJ^ii—^yS- in Hq, and subtract it from V, and adjust the strength t to get best convergence 
in the series. The basic idea of the method has been discussed before so we only give 
brief details here. To compute series for C(r) in powers of A we add a field term to H 

H = Ho + XV + hY,S,-S,+r (6) 

i 

compute the ground state energy in the form 

£oiX,h) = EoiX) + hNCir)/A + Oih'') (7) 

and hence extract series in A for C{r). For the longitudinal correlator the field term is 
hY,i SfSfj^^. Examples are give in the following sections. 
Two important questions are not addressed in this work: 

(i) the behaviour of correlators at large distances, and the asymptotic behaviour and 
associated critical exponents. 

(ii) correlators at finite temperature. 

II. THE 1-DIMENSIONAL CASE 

The 5* = I antiferromagnetic Heisenberg chain is exactly solvable by the Bethe ansatz 
0,^, and the ground state energy and elementary excitations are given by simple analytic 
expressions. However the wavefunction is sufficiently complex that little exact information 
is available on correlators. In fact only the first two are known exactly, and are 

(Si ■ S,+i) = ^(1 - 41n2) = -0.443147... 

(8) 

(Si ■ Si+2) = ^(1 - 161n2 + 9C(3)) = 0.182039... 

The first result comes from the ground state energy , while the second is obtained via 
the strong coupling limit of the Hubbard model. There is no spontaneous symmetry break- 

3 



ing in the ground state of the isotropic spin chain, and hence longitudinal and transverse 
correlators are equal, and obtainable directly from (|^). 

The first serious attempts to obtain further results for the antiferromagnetic chain were 
by Bonner and Fisher [0, who used exact diagonalizations for systems up to = 10 spins, 
and by Kaplan and co-workers |Q who extended this to = 18. These results suffer from 
large finite-size effects, and need to be extrapolated to the therodynamic limit via a finite- 
size scaling ansatz. In this way Kaplan et al. estimated the value of correlators up to 8th 
neighbors, with confidence limits of about 1% in C(8). We show these values in Table |. 
Subsequently Lin and Campell extended the exact diagonalizations to = 30. By use 
of the empirical scaling relation 

CN{r) = CUr)f{r/N) (9) 

with 

f{y) = [1 + 0.28822 sinh^ (1.673?/)] ^'^^ (10) 

they estimated correlators up to r = 15, i.e. 15th neighbors. However the accuracy of this 
scaling is perhaps doubtful since it is known that there are logarithmic terms which slow 
convergence |]10| . 



The development of the Density Matrix Renormalization Group (DMRG) method |TT 



allows much longer chains to be treated with high numerical accuracy. Hallberg et al. |T2 
have used DMRG to compute correlators for Heisenberg chains up to = 70 spins, with 
a scaling function similar to (|^) used to extrapolate to the thermodynamic limit. The data 
were shown to be consistent with the asymptotic behaviour 

C{r) ^ {-lyihirf'^/r (11) 



predicted by field theory pO |. 

We have described the series method briefly in the Introduction. Using this approach we 
have computed expansions in A, for both the total and longitudinal correlators for distances 



r = 1, 2, ■ ■ ■ , 10. The maximum order is A^'* for r = 1 and A^^ for r = 10. We note that 
the longitudinal correlators, and the total correlators for r even, contain only even powers 
of A. We also note that the series are rather erratic, both in sign and magnitude of the 
coefficients. This had already been noted by Walker [|13| who expanded the ground state 
energy, and hence C(l), to order A^^. Rather than quote all series here we will make them 
available to any reader on request. Table |T| shows the coefficients for the series for C(4), 
C/(4) and Cf(4). We note that, as expected, the series for the transverse correlator Ct{r) 
starts with a term A*". 

The series have been evaluated for fixed A by means of integrated differential approx- 
imants [1^, and the values of correlators for r = 1,2, ■■■,6 are shown in Figure 1. The 



analysis becomes less precise as the weakly singular point A = 1 is approached. We also 
show in the Figure the extrapolated exact diagonalization results. As can be seen from the 
Figure, and from Table ||, the agreement is very good. It is clear that in 1-dimension the 
series method is not able to match the precision of either scaled finite lattice or DMRG 
results, but in higher dimension these latter methods are not competitive. Furthermore, as 
we shall show, the series analysis can be made more precise in d > 2 because the stronger 
singularity at A = 1 can be removed by a transformation, and the Ising expansions used 
here are more suitable for d >2, where the ground state has long range Neel order. 



We should mention here also the work of Singh et al. [|T5| who used exactly the same 
method as ours to compute the structure factors 

-| oo 

S.. = - E [{-irC,{r)~A<Sl>'] (12) 

r=— oo 

and 

-| oo 

-5+- = 2 E i-^rC,{r) (13) 

r=— oo 

for the S = \ antiferromagnetic chain. Our correlator series, when summed, agree with their 
results. 
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III. THE 2-DIMENSIONAL CASE 



There has been much interest, in recent years, in the nature of the ground state of the 
Heisenberg antiferromagnet on the square lattice. There is considerable evidence, from exact 
diagonalizations |ll6|-[T^ and quantum Monte Carlo calculations [|1^-|21]] that the ground state 



breaks rotational symmetry, giving rise to a staggered magnetization in some direction. This 
is generally referred to as a quantum Neel state, with Neel type order reduced to approx. 
60% of its classical value by quantum fluctuations. The situation is summarized in recent 
reviews 



p2| , p3| . The focus has generally been on the ground state energy and staggered 
magnetization, although some short range correlators have also been computed. 

In any finite system there can be no spontaneous symmetry breaking and hence the 
exact diagonalization and Monte Carlo studies cannot distinguish between the longitudinal 
and transverse correlators for the isotropic case. Furthermore if Ci or Ct are computed by 
these methods the values will not yield correct results for the thermodynamic limit, where 



Other approaches, such as spin wave (SW) theory [p^ -|27[], variational methods ||28|, or 
perturbation series about the Ising limit [^^^ start from a broken symmetry state, which 
is preserved during the calculation. It seems highly likely, although we know of no proof, 
that these approaches will yield the correct symmetry-broken state of the infinite isotropic 
system. 

We have computed series expansions for a number of local correlators for the square 
lattice S = ^ antiferromagnet. The expansions start from the Ising limit and are carried 
through order 14,9,9,9,7 for C(r), C/(r) with r = (0, 1), (1,1), (2,0), (3,0), (4,0) The series 
coefficients (for t = 0) are given Table |ITT[ 

In analysing the series it is advantageous to transform to a new variable 

5 = 1 - (1 - A)i/2 ^ (14) 

to remove the singularity at A = 1. Spin wave theory predicts a square root singularity 
of this type. This transformation was first proposed by Huse |^ and was also used in 



earlier work on the square lattice case p^. We then use both integrated first-order inho- 



mogeneous differential approximants [0] and Pade approximants to extrapolate the series 
to the isotropic point 5 = 1 (A = 1). The results are shown as functions of A in Figure 2 for 
r=(l,0), (1,1). We also show the transverse correlator, obtained from Eq. In the Ising 
limit the total and longitudinal correlators are equal and the transverse correlator is zero. 
As we increase the transverse coupling, the longitudinal correlators decrease in magnitude 
while, as expected, the transverse correlators increase, while the total correlator increases 
in magnitude for nearest neighbors, but is reduced for second neighbors. The behaviour of 
further correlators is similar, and is not shown. It is also noteworthy that at A = 1 the lon- 
gitudinal and transverse correlators remain unequal, reflecting the symmetry broken ground 
state. In Fig. 3 we show a comparison between our series results and other methods for 
the nearest neighbor correlators. For small A all methods are in close agreement, but near 
the isotropic point linear spin wave theory become poor for longitudinal (and transverse) 
correlators, whereas exact finite lattice diagonalizations have longitudinal and transverse 
correlators equal at A = 1. Third order spin- wave theory is much better, being almost 
indistinguishable from the series results over the whole range of A. 

In Table we give numerical estimates of all the correlators at the isotropic point, 
obtained by our series method and by exact diagonalization/Monte Carlo on finite lattices 



17| , |20[| and linear spin wave theory 0. We believe that the expressions in Ref. [^] contain 



minor errors, and should read, for and r on the same sublattice 



(15) 



^ k ^1 - XH 

while for and r on different sublattices: 



(16) 



where the notation is as in Ref. and A is the anisotropy parameter. 



We note from Eqn. |1^ and |1^ and Table |^ that first order spin-wave theory gives a 
longitudinal correlator which is independent of distance, clearly an artifact of the approxi- 
mation. The total correlator is however very consistent with the series results. The picture 
is considerably improved in higher order spin-wave theory where, for example, 3rd order 
spin wave theory gives 3-figure agreement with series for all of C, Ci and Ct for nearest 
neighbors. We have not attempted to carry this out for further neighbors, and are unaware 
of any work along these lines. 

IV. THE 3-DIMENSIONAL CASE 

We have used the same series approach to calculate correlators for the S = ^ antifer- 
romagnet on the simple cubic lattice. The magnetically ordered ground state will again 
be reflected in a difference between longitudinal and transverse correlators at the isotropic 
limit. 

Expansions, starting from the Ising limit, have been obtained for C(r), Ci{r) for the five 
cases r =(1,0,0), (1,1,0), (2,0,0), (3,0,0), (4,0,0) to order 12,7,9,7,7 respectively. We have 
again used a staggered field term t J2i{~^yS^ to improve convergence. The series coefficients 
(for t = 0) are given in Table 0. The series is extrapolated in a similar way as that for the 
square lattice. Figure 4 shows the nearest and next-nearest neighbor correlators as functions 
of the anisotropy parameters. This is qualitatively similar to Figure 2, but clearly shows that 
in 3-dimensions transverse correlators are reduced and the difference between transverse and 
longitudinal correlators is increased for all values of the anisotropy parameter. 



In Table ^ we give numerical estimates of all correlators at the isotropic point, and 
a comparison with 1st order spin-wave theory. It is apparent that the correlators fall off 
more slowly with distance than in the 2-dimensional case, reflecting the greater stability of 
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antiferromagnetic long-range order in the ground state in 3-dimensions. It is also apparent 
that the transverse correlators are, relatively, much weaker in 3-dimensions, consistent with 
weaker quantum fluctuations. Linear spin-wave theory gives reasonable results for the total 
correlators, but again suffers from the defect of having longitudinal correlators independent 
of distance. Third-order spin-wave theory gives results for nearest neighbor correlators in 
excellent agreement with the series results. 



V. CONCLUSIONS 

We have used series methods to obtain numerical estimates for short-distance ground 
state correlation functions for the S = j Heisenberg antiferromagnet on square and simple 
cubic lattices. Despite their importance in characterising the nature of the antiferromagnetic 
ground state, there appears to have been little previous work on the subject. 

The series approach is able to provide rather precise estimates for correlators up to at 
least 4 lattice spacings. The results reflect the known breaking of rotational symmetry in 
the ground state, in that longitudinal and transverse correlators remain unequal even in 
the isotropic Hamiltonian limit. Exact diagonalizations and Monte Carlo calculations on 
finite lattices are unable to account for this and hence will not yield correct estimates for 
longitudinal and transverse correlators separately. In 3-dimensions no results are available 
from diagonalizations or quantum Monte Carlo, beyond nearest neighbors. 

We have shown that first-order spin wave theory gives rather poor estimates but 2nd and 
3 rd order spin wave theory gives excellent agreement with series results for nearest-neighbor 
correlators. Higher order spin wave results have not been obtained for further correlators, 
to our knowledge. 

As a test of the method we also computed correlation series for the 1-dimensional case. 
The results were quantitatively accurate, but less precise than the DMRG method. 

This approach can also be used to calculate correlators for more complex models involving 
competing interactions. For example, we have studied the Ji — J2 model PT|, which has a 
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quantum critical point at J2/J1 ^ 0.38, where the Neel order is destroyed and the system 
enters a magnetically disordered spin-liquid phase. We find that the difference between 
longitudinal and transverse correlators remains nonzero in the Neel phase, but vanishes at 
the quantum critical point, indicating a restoration of full rotational symmetry in the ground 
state at that point. We expect this method to prove useful in other problems of this type. 
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TABLES 



TABLE L Estimates of the total correlators C(r) for the isotropic -S* = 2 antiferromagnetic 
chain (Eq. ^ in text) from exact diagonalizations Q and series (the present work). 



r 


Exact Diagonalization 
+ rinite oize ocanng (^nei. |o||j 


Series Expansion 
evaluated at A = 1 


1 


1 770/1 /'f;^ 


-i. M 6\Z ) 





U. (90(^0 j 


n 7Q^/'Q^ 
U. / oU(o j 


3 


-0.6027(3) 


-0.588(15) 


4 


0.4158(10) 


0.408(20) 


5 


-0.3705(10) 


-0.38(3) 


6 


0.2946(12) 


0.32(4) 


7 


-0.2697(20) 




8 


0.2280(20) 
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TABLE II. Coefficients of correlator series for the S = ^ antiferromaj 



gnetic chain for r = 4 



power of A 



C(4) 



a(4) 




2 
4 



10 
12 
14 
16 
18 
20 
22 



1.000000000 
-2.000000000 

2.000000000 
2.500000000x10"^ 

-1.265625000 
3.281250000x10"^ 

-4.296875070 xl0"2 
1.225585910x10"^ 
6.085209298x10"^ 
5.311830238xl0"3 

-9.656318929x10-3 



1.000000000 
-2.000000000 

5.000000000x10"^ 
2.375000000 

-2.250000000 
5.546875000x10"^ 

-1.865234375 xlO"i 
2.441406766x10-2 
7.099156584x10"^ 
4.902810212x10-2 
2.493706992 xl0"2 
7.884341059xl0"3 



0.000000000 
0.000000000 

7.500000000x10"^ 
-1.062500000 

4.921875000x10"^ 
-1.132812500x10-1 

7. 177734342 xl0"2 

4.907226169x10-2 
-5.069736427x10-3 
-2.185813594x10-2 
-1.729669442x10-2 
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TABLE III. Non-zero coefficients of various correlator series for the S = ^ antiferromagnet on 



the square lattice {t = 0) 



power of A 


C(r) 


Ci{r) 


Ct{r) 





-1.000000000 


r=(l,0) 

-1.000000000 


0.000000000 


1 


-6.666666667x10-1 


0.000000000 


-3.333333333x10-1 


2 


3.333333333x10-1 


3.333333333x10-1 


0.000000000 


3 


7.407407407x10-3 


0.000000000 


3.703703704x10-3 


4 


-5.555555556x10-3 


-5.555555556x10-3 


0.000000000 


5 


-1.897883598x10-2 


0.000000000 


-9.489417989x10-3 


6 


1.581569665x10-2 


1.581569665x10-2 


0.000000000 


7 


-1.320340554x10-2 


0.000000000 


-6.601702770x10-3 


8 


1.155297985x10-2 


1.155297985x10-2 


0.000000000 


9 


-6.237012985x10-3 


0.000000000 


-3.118506492x10-3 


10 


5.613311689x10-3 


5.613311689x10-3 


0.000000000 


11 


-5.806609913x10-3 


0.000000000 


-2.903304957x10-3 


12 


5.322725757x10-3 


5.322725757x10-3 


0.000000000 


13 


-4.231435003x10-3 


0.000000000 


-2.115717502x10-3 


14 


3.929189659x10-3 


3.929189659x10-3 


0.000000000 





1.000000000 


r=(l,l) 

1.000000000 


0.000000000 


2 


-2.222222222x10-1 


-4.444444444x10-1 


1.111111111x10-1 


4 


3.444444444x10-2 


4.567901235x10-3 


1.493827160x10-2 


6 


7.314352902x10-"' 


-1.979035301x10-2 


1.026089415x10-2 


8 


-1.582060575x10-3 


-1.475015802x10-2 


6.584048724x10-3 





1.000000000 


r=(2,0) 

1.000000000 


0.000000000 


2 


-3.333333333x10-1 


-4.444444444x10-1 


5.555555556x10-2 


4 


4.246913580x10-2 


-1.666666667x10-2 


2.956790123x10-2 


6 


3.666832535x10-3 


-1.353930181x10-2 


8.603067173x10-3 


8 


-3.706741518x10-3 


-1.786067634x10-2 


7.076967412x10-3 





-1.000000000 


r=(3,0) 

-1.000000000 


0.000000000 


1 


0.000000000 


0.000000000 


0.000000000 


2 


4.444444444x10-1 


4.444444444x10-1 


0.000000000 


3 


-3.888888889x10-2 


0.000000000 


-1.944444444x10-2 


4 


2.037037037x10-2 


2.037037037x10-2 


0.000000000 


5 


-4.035089653x10-2 


0.000000000 


-2.017544827x10-2 


6 


2.054058327x10-2 


2.054058327x10-2 


0.000000000 


7 


-1.401574252x10-2 


0.000000000 


-7.007871260x10-3 
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8 1.785822451x10-2 1.785822451x10-2 0.000000000 

r = (4, 0) 

1.000000000 1.000000000 0.000000000 

2 -4.444444444x10-1 -4.444444444 X 10" ^ 0.000000000 

4 -1.172839506x10-2 -2.172839506x10-2 5.000000000 x 10-^ 

6 -1.337086028x10-3 -2.198640296x10-2 1.032465846x10-2 



TABLE IV. Values of correlators for the isotropic S = ^ Heisenberg antiferromagnet on the 



square lattice. 



r 




Series (this work) 




Finite Lattice^ 


Linear Spin Wave Theory'' 




C{r) 


Ci{r) 


Ct{r) 


C(r) 


C(r) 


Ci{r) 


Ct{r) 


(1,0) 


-1.3386(2) 


-0.572(4) 


-0.383 


-1.344, - 


-1.316 


-0.2136 


-0.551 


(1,1) 


0.794(15) 


0.430(6) 


0.182 


0.765, 0.82 


0.795 


0.2136 


0.291 


(2,0) 


0.67(2) 


0.408(10) 


0.131 


0.84, 0.71 


0.673 


0.2136 


0.230 


(3,0) 


-0.52(2) 


-0.386(10) 


-0.067 


-0.75, -0.60 


-0.526 


-0.2136 


-0.156 


(4,0) 


0.42(2) 


0.376(20) 


0.022 


-, 0.53 


0.440 


0.2136 


0.113 


^Prom exact diagonalizations, N = 


26 (Ref. 


flTll), and project' 


or Monte Carlo (Ref. Q). 




^For r = 


(1,0) we 


also have result 


s from 2nd and 3rd order 


spin- wave theory (Ref. [p7|) 


, which 


give for C 


, C/ and Ct respectively - 


-1.3408,- 


-0.672, -0.334 (2nd order) and 


-1.3400, -0.575 


, -0.383 



(3rd order). 
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TABLE V. Non-zero coefficients of various correlator series for tfie S = ^ antiferromagnet on 



tfie simple cubic lattice {t = 0) 



power of A 


C(r) 


Ci{r) 


Ct{r) 





-1.000000000 


r = (1,0,0) 

-1.000000000 


0.000000000 


1 


-4.000000000x10-1 


0.000000000 


-2.000000000x10-1 


2 


2.000000000x10-1 


2.000000000x10-1 


0.000000000 


3 


2.666666667x10-3 


0.000000000 


1.333333333x10-3 


4 


-2.000000000x10-3 


-2.000000000x10-3 


0.000000000 


5 


-1.252256756x10-2 


0.000000000 


-6.261283778x10-3 


6 


1.043547296x10-2 


1.043547296x10-2 


0.000000000 


7 


-5.018108039x10-3 


0.000000000 


-2.509054019x10-3 


8 


4.390844534x10-3 


4.390844534x10-3 


0.000000000 


9 


-3.490020757x10-3 


0.000000000 


-1.745010379x10-3 


10 


3.141018683x10-3 


3.141018683x10-3 


0.000000000 


11 


-2.400011436x10-3 


0.000000000 


-1.200005718x10-3 


12 


2.200010484x10-3 


2.200010484x10-3 


0.000000000 





1.000000000 


r = (1,1,0) 

1.000000000 


0.000000000 


2 


-1.600000000x10-1 


-2.400000000x10-1 


4.000000000x10-2 


4 


1.766543210x10-2 


7.567901234x10-" 


8.454320988x10-3 


6 


-3.636588526x10-3 


-1.226438856x10-2 


4.313900018x10-3 





1.000000000 


r = (2,0,0) 

1.000000000 


0.000000000 


2 


-2.000000000x10-1 


-2.400000000x10-1 


2.000000000x10-2 


4 


1.251358025x10-2 


-1.106172840x10-3 


6.809876543x10-3 


6 


-4.632761555x10-3 


-1.262448646x10-2 


3.995862454x10-3 


8 


-9.042568955x10-* 


-5.558061905x10-3 


2.326902505x10-3 





-1.000000000 


r = (3,0,0) 

-1.000000000 


0.000000000 


1 


0.000000000 


0.000000000 


0.000000000 


2 


2.400000000x10-1 


2.400000000x10-1 


0.000000000 


3 


-8.148148148x10-3 


0.000000000 


-4.074074074x10-3 


4 


1.550617284x10-3 


1.550617284x10-3 


0.000000000 


5 


-5.458926876x10-3 


0.000000000 


-2.729463438x10-3 


6 


1.292980049x10-2 


1.292980049x10-2 


0.000000000 


7 


-4.261178891x10-3 


0.000000000 


-2.130589445x10-3 





1.000000000 


r = (4,0,0) 

1.000000000 


0.000000000 


2 


-2.400000000x10-1 


-2.400000000x10-1 


0.000000000 


4 


-4.098765432 xlO-"! 


-1.644444444x10-3 


6.172839506x10-"* 
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6 



-1.131207676x10-2 



-1.297046176x10-2 



8.291924985x10-4 



TABLE VI. Values of correlators for the isotropic S = 2 Heisenberg antiferromagnet on the 



simple cubic lattice. 



r 




Series (this work) 






Linear Spin Wave Theory'' 




C(r) 


C'lir) 


Ctir) 


C(r) 


Ci{r) 


Ct{r) 


(1,0,0) 


-1.2028(3) 


-0.775(3) 


-0.214 


-1.1943 


-0.6866 


-0.2539 


(1,1,0) 


0.857(2) 


0.683(8) 


0.087 


0.8440 


0.6866 


0.0787 


(2,0,0) 


0.807(2) 


0.684(8) 


0.061 


0.7900 


0.6866 


0.0517 


(3,0,0) 


-0.768(8) 


-0.676(8) 


-0.046 


-0.7317 


-0.6866 


-0.0226 


(4,0,0) 


0.755(8) 


0.672(8) 


0.041 


0.7097 


0.6866 


0.0116 



^For r = (1,0,0) we also have results from 2nd and 3rd order spin-wave theory (Ref. [^), which 
give for C, C/ and Q respectively -1.2038,-0.7756,-0.2141 (2nd order) and -1.2033, -0.770(2), 
-0.2165(9) (3rd order). 
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FIGURES 

FIG. 1. Correlators C(r) for the S = ^ antiferromagnetic chain for r = 1,2, -6. The full lines 
give the total correlator C(r), the dashed lines give 3Ci{r). Curves for different r are labelled at 
the right hand edge. The circles at A = 1 are the values from the finite-lattice calculations 
Note that at the isotropic point A = 1 the longitudinal and transverse correlator are equal. 

FIG. 2. Correlators for nearest and next-nearest neighbors for the square lattice, for varying 
anisotropy parameter A. Full lines denote the total correlator C(r), dashed lines the longitudinal 
correlator C/(r), and dotted lines the transverse correlator Ct(r). 

FIG. 3. Comparison between series (this work) and other estimates of the nearest neighbor 
correlator for the square lattice, for varying anisotropy parameter A. 

FIG. 4. Correlators for nearest and next-nearest neighbors for the simple cubic lattice. Full, 
dashed and dotted lines represent total, longitudinal and transverse correlators respectively. 
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